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Keywords Abstract. This paper investigates different liquidity management 
Islamic Banking aspects of Islamic banking industry in Pakistan. Islamic banking liabilities 
Financial Intermediation (deposits) and assets models, with regard to liquidity management, have 
Risk Management found the significant role of following variables: (a) returns on deposits, (b) 
Liquidity Management returns on financing, (c) costs of banking operations, and (d) the Interbank 
Auto Regressive Distributed rate, here Karachi Interbank Offer Rate (KIBOR). Islamic banking liquidity 
Lag (ARDL) Modelling reserves model, however, recommends that Islamic banks need to consider 
ARIMA Modelling following variables, while developing optimum liquidity reserves: (a) total 


Islamic financing, (b) returns on financing, and (c) KIBOR. Moreover, 
the resilience analysis of the Islamic banking industry carried out in the 
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basis of Autoregressive Integrated Moving Average (ARIMA) models for 
the current study suggest that tier-2 liquid instruments would possibly be 
performing well in mitigating any future liquidity run conditions (up to 
95% of deposits). However, in case of tier-1 liquid instruments, there is 
a possibility of liquidity mismatches when liquidity withdrawals exceed 
the limit of 55% of deposits. In conclusion, Islamic banking depositors, 
besides their religious motives of supporting Islamic banks, expect from 
their banks to earn profits and pay competitive returns on their deposits. 
Therefore, Islamic banks need to make prudent portfolio financing so as to 
pay competitive returns to their depositors. 
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INTRODUCTION 


Liquidity management means strengthening a bank’s ability to meet all its financial obliga- 
tions on its liability side, as well as, seizing all the investment opportunities on its assets’ 
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side, without incurring any unexpected losses (Basel Committee, 2008). According to Ac- 
counting and Auditing Organization for Islamic Financial Institutions (AAOIFI’s) Shari‘ah 
standard on liquidity management, liquidity refers to the money and other financial assets 
that are easily convertible into money; whereas, liquidity management, in case of a bank, 
refers to the realization of an appropriate and befitting balance between liquidity’s acqui- 
sition, at the best price in the shortest possible time, and liquidity’s investment and/or em- 
ployment in the best possible manner (AAOIFI, 2014). 

A bank requires liquidity to accommodate every fluctuation, expected as well as unex- 
pected, of its balance sheet. That is, a bank is required to have ample liquid funds: (a) to 
fulfil depositors’ withdrawal demands without any delay; (b) to pay maturing liabilities; (c) 
to fulfil financing demands by entrepreneurs; and (d) to rebalance the investment portfolio 
(Saiti, Hasan, & Engku-Ali, 2016). A bank will be experiencing the liquidity risk, which is 
the spectre haunting the banking industry, when it would not receive liquid funds in times of 
dire needs, except with difficulty and at high costs (Wahyudi, Rosmanita, Prasetyo, & Putri, 
2015). 

Despite impressive development of Islamic finance over last three decades, a complete 
Islamic financial system, with its distinctive financial instruments and financial markets, 
is still in early stages of evolution. It is pertinent to note that Islamic banking system is 
comparatively well-developed; while, Islamic capital and money markets are at their infancy 
stages and experiencing a shortage of Shari‘ah compliant tools used in both mobilization 
and utilization of funds (Ahmed, 2011; Zaher & Hassan, 2001). According to some recent 
estimates, total assets held by the Islamic financial institutions worldwide are about USD 
1.8 trillion (Rizvi & Alam, 2016). Whereas, the total assets held by Pakistan-based Islamic 
financial institutions have reached to about PKR 1.7 trillion (SECP, 2015). 

Roots of Islamic finance can be traced back to the early Islamic history and beyond. It 
is reported that trade and commerce in the Medieval Islam were carried out on the basis 
of partnership (musharakah), commenda (muddrabah), credit sale (bai‘ al-mw ajjal), ad- 
vance payment sale (bai‘ al-salam), and benevolent lending (gard al-hasan) (Goitein, 1971; 
Udovitch, 1970). It is worth mentioning here that Islamic financial system is prohibition- 
driven. According to Ibn-e-Taymiyya (d. 728 A.H./1328 A.D.), two prohibitions, that of 
riba (interest) and gharar (uncertainty with regard to the subject matter and the price in 
contracts), can distinguish between financial contracts that are valid or void (El-Gamal, 
2006). The aforementioned prohibitions, along with the prohibition of maysir (gambling- 
like speculations) and regard for others’ interests while pursuing material gains, are central 
to the Islamic economic teachings (Siddiqi, 2006; 2009). 

It is advocated by Islamic economists that: since, Islamic finance is based on real eco- 
nomic activities, and devoid of excessive leverages and imprudent risk taking; therefore, it 
has the ability to correct the structural flaws of the conventional finance due to which finan- 
cial crises have been occurring time and again (SBP, 2015). Quaid-e-Azam Muhammad Ali 
Jinnah, Muslims can fulfil their mission of giving humanity a message of peace only by the 
development of banking practices in the light of Islamic principle of socio-economic justice 
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he rejected the western economic theory and finance system and emphasized Muslims to 
present an Islamic economic system promoting equality and socio-economic justice’. 

According to financial intermediation theory, the maturity transformation of short-term, 
liquid deposits into long-term, illiquid loans (financing in case of Islamic banks) is the 
fundamental role of banks, which inherently makes banks vulnerable to liquidity risk (Basel 
Committee, 2008; Berger & Bouwman, 2009). The primary aim of the research paper is to 
analyse econometrically different liquidity management aspects of Islamic banking industry 
in Pakistan. The aim is further broken down into the following two objectives: (i) to identify 
factors that influence Islamic banks in managing liquidity on liability as well as asset sides 
and in maintaining the optimum liquidity reserves; and (11) to investigate the resilience of 
the Pakistan-based International Banking Facility (IBF) industry against certain liquidity 
pressures. 

The rest of the paper is organised as follows: Section 2 provides a brief literature re- 
view on liquidity management by banking institutions, while section 3 explains the research 
methodology of the current study. Afterwards, section 4 analyses econometrically differ- 
ent liquidity management aspects of Pakistan-based IBF industry and develops econometric 
models, and finally, section 5 concludes the paper. 


LIQUIDITY MANAGEMENT BY ISLAMIC BANKS 


Islamic finance has a long history going back to the early days of Islam and its theories are 
derived from revealed texts. Islam encourages its believers to conduct trade and business 
being fair, honest and just towards others. It is pertinent to note that prohibitions of riba 
and gharar most importantly differentiate Islamic business and commerce from western 
commercial transactions (Venardos, 2012). 

It has been reported that the Prophet (S.A.W) himself participated in mudarabah con- 
tracts as mudarib, while capital providers (rab al-m@l) were his wife and others. Moreover, 
his companions of such as Umar (R.A) and Usman (R.A), used to invest funds belonging to 
orphans on the basis of mudarabah. It is also reported that Muslim traders were using finan- 
cial instruments-sukak-in order to mitigate risks involved (theft or loss of money) in carry- 
ing money on long commercial journeys between Islamic world and East Asia (Jamaldeen, 
2012). 

Financial intermediation in the medieval Islamic world was closely tied to commercial 
activities. It is pertinent to note that sarraf (goldsmiths and money exchangers) of that 
time were supporting commercial activities by undertaking various traditional functions 
of a modern financial institution, such as, operating secure cross-border payment systems 
through sukuk issuance (Udovitch, 1981). Interestingly, sarraf were operating through or- 
ganized networks that established them as intermediaries of the medieval period (Chapra & 
Ahmed, 2002). It has also been reported that some of the contracts, concepts and institutions 
developed in the medieval Islamic world of late eighth century provided foundations for the 





‘For details, see the inaugural address of Quaid-e-Azam Muhammad Ali Jinnah available at: 
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development of similar contracts and instruments in the Western world several centuries 
later (Udovitch, 1970). 

Islamic banks typically attract funds in the form of: (a) Current Accounts deposit holders; 
and (b) unrestricted Investment Account Holders (IAHs). It is pertinent to note that C/As 
and some demand deposits are in the form of gard al-hasan or wadi‘ah to Islamic banks, 
which can be withdrawn anytime. For the purpose, Islamic banks need to maintain a sound 
repayment capacity to meet fully the cash withdrawals by demand deposit holders. Unre- 
stricted investment accounts, on the other hand, are investments made by IAHs on the basis 
of bai‘ al-mw’ajjal. Therefore, IAHs share in profits and losses (Al-Amine, 2013; Central 
Bank of Oman, 2012; Iqbal & Mirakhor, 2011; Vogel & Hayes, 1998). 

Liquidity management is a shared responsibility of Islamic banks and their regulators, 
i.e., central banks. Islamic banks, like conventional banks, generally maintain their liquidity 
positions to such an extent that facilitate them in: (a) meeting their commitments; (b) main- 
taining required reserves; and (c) seizing all the investment opportunities-while remaining 
solvent and profitable. It is important here to realize that liquidity and solvency-in case of 
banks (whether Islamic or Conventional)-are interlinked: an illiquid bank may soon become 
insolvent, while, an insolvent bank is illiquid, and therefore, causes the financial instabil- 
ity. Central banks, on the other hand, are equally responsible for liquidity management of 
banks, in order to: (a) ensure an efficiently stable banking system; (b) ensure the smooth 
interbank payment settlements; and (c) maintain a level of liquidity in the banking system 
that is consistent with the central bank’s monetary policy targets (Ganley, 2004; Kahf & 
Hamadi, 2014). 

Islamic banks worldwide rely heavily on sukik-an Islamic alternative to conventional 
bonds-for the liquidity management purposes. It is pertinent to note that the demand for 
sukuk is much greater than their availability. This demand-supply gap of sukik has left no 
other choice for Islamic bankers to use Commodity Murabahah-as a liquidity management 
instrument. Commodity Murabahah (based on tawarrugq) is controversial in terms of its 
Shart‘ah permissibility. Some contemporary Shart‘ah scholars have opined that Commod- 
ity Murabahah should only be used in extreme emergency only to avoid dealing in interest 
(Usmani, 2015). However, Kahf and Hamadi (2014) argued that the commodity Murabahah 
is not permissible since this transaction involves an organized banking tawarruq, which is 
not permitted in the Shart‘ah perspective as resolved by the Organisation of Islamic Coop- 
eration (OIC) Figh Academy. 

Pakistan, with the world’s 2™ largest Muslim population, is the first Muslim country 
where, during 1950s, an Islamic mutual savings bank was established by religious-minded 
landlords to provide agricultural finance to poor farmers (Asutay, Aysan, & Karahan, 2013; 
Siddiqi, 2006), while, the Islamisation process of the whole economy was initiated in late 
1970s. But such a challenging task still remains incomplete due to lack of interest of dif- 
ferent stakeholders, apart from some other reasons. In the meanwhile, Iran and Sudan got 
established their interest-free banking systems. As per the amended policy, in Pakistan, the 
first Islamic bank was given the license and became operational in 2002 (Beseiso, 2014; 
Khan & Bhatti, 2008; Farooq & Zaheer, 2015). 


2018 Journal of Islamic Business and Management Vol. 8 Issue 1 57 


Presently, there exists a dual banking system in Pakistan, where conventional and Islamic 
banks coexist and compete with each other. It is worth mentioning here that, during the fi- 
nancial panic of 2008, Pakistan-based Islamic banks experienced fewer deposit withdrawals 
as compared to their conventional counterparts; even some Islamic banks recorded deposit 
increases, which lead to a net deposits’ inflow into Islamic banks. Thus, Islamic banks have 
shown superior performance, greater resilience, and at the same time, ability to inject liq- 
uidity into the real economy during times of financial distress as well (Farooq & Zaheer, 
2015). 

In a recent study titled, “Knowledge, Attitude, and Practices of Islamic banking in Pak- 
istan” (KAP-study), jointly conducted by State Bank of Pakistan (SBP) and Department 
for International Development (DFID), U. K., an overwhelming demand for Islamic bank- 
ing has been identified in the Pakistan. The findings of the KAP-study further reveal that 
Islamic banking demand is higher amongst households (95%) as compared to businesses 
(73%). Moreover, 88.41% Banked and 93% Non-Banked respondents of the KAP-study 
considered that the interest charged and given by banks is the prohibited riba (SBP, 2014). 

After the re-launch since 2002 Islamic banking, with the full support of the government, 
has now become the fastest growing segment of the Pakistan-based banking industry. There 
is a widespread agreement among the academicians and practitioners of Islamic banking 
and finance in Pakistan that finding an innovative, Shart‘ah compliant solution to the liquid- 
ity management challenge of Islamic banks, together with increasing the awareness in the 
country regarding the Islamic banking model, are the two important tasks ahead (Russell- 
Walling, 2014). 

In Pakistan, Islamic banks are required to maintain the same Cash Reserve Requirement 
(CRR). However, the Statutory Liquidity Requirement (SLR) for Islamic banks has been set 
differet from that of conventional banks since 2017. Conventional banks have opportuni- 
ties of investing in Pakistan Investment Bonds (PIBs) as well as treasury bills (T-bills). Is- 
lamic banks on the other hand, are lacking a Shart‘ah-compliant alternative of T-bills (Alam, 
2014). Under such circumstances, Pakistan-based Islamic banks rely heavily on sukuk for 
investing their excess liquidity. Sukik are also held by Islamic banks for SLR purposes. 
According to some recent estimates, total sukik issued in Pakistan amount to nearly PKR 
695 Billion (Saheed, 2014), which are not enough to meet the Islamic banks’ demand. 

Islamic banks in Pakistan generally use two interbank musharakah and interbank wakalah 
for short-term liquidity placements. Some Islamic banks are also using Commodity Mura 
bahah for investing their excess liquidity. 


RESEARCH METHODOLOGY 


The research methodology of this paper is quantitative. The paper has developed econo- 
metric models that explain different liquidity management aspects of Pakistan-based IBF 
industry. For this purpose, secondary data has been collected from State Bank’s website and 
Islamic banks’ financial statements. 

Quantitative research explains phenomena by collecting numerical data associated with 
the phenomena and analyse it using statistical methods (Aliaga & Gunderson, 2002). It 
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typically involves: (i) formation of hypothesis from theory; (ii) collection of primary or 
secondary data associated with the social phenomenon under investigation; (iii) analysing 
the numerical data either statistically or econometrically in order to test the hypothesis; (iv) 
producing the research’s output explaining the social phenomenon under investigation (Is- 
mal, 2010). 


ARDL Models 

An ARDL model regresses a dependent variable over the present and past values of a set 
of independent variables along with the past values of the dependent variable (Fabozzi, 
Focardi, & Kolm, 2006). This paper uses ARDL approach to model Islamic banks’ bal- 
ance sheet, because of the following: (i) each variable of Islamic banks’ balance sheet can 
possibly be functioning in an equation either as an independent variable or as a dependent 
variable; (ii) other than present values of independent variables, their time lags (past values) 
are also explaining the dependent variable in an ARDL model; and (iii) an ARDL model 
provides a vehicle for testing presence of long-run relationships among variables (Ismal, 
2010) 

ARDL models, estimated using Ordinary Least Squares (OLS), are most commonly used 
econometric tools for developing robust dynamic models. Other econometric tools, for 
the purpose of estimating robust dynamic models, require an extensive series of data that 
is not available in case of Islamic banks in Pakistan. It is mainly due to the reason that 
Pakistan-based Islamic banking industry is still in an evolutionary stage with limited data 
and information. 

Time series data of Islamic banking variables, ranging between July 2006 and December 
2016, is obtained from SBP’s website and quarterly financial statements of Islamic banks in 
Pakistan. It is pertinent to note that the modelling process follows the following sequential 
steps in order to make robust dynamic models: (i) developing a theoretical model on the 
basis of literature review; (ii) selecting independent variables; (iii) selecting hypothetically 
signs of slope coefficients; (iv) collecting, inspecting and cleaning the data; (v) estimating 
the regression equation and evaluating it; and (vi) documenting regression results. (Student- 
mund, 2016). 


ARIMA Models 
ARIMA modelling techniques are being used in the current study to: (a) estimate models 
for variables related to Islamic bank’s liquidity management by employing their historical 
data from July 2006 to December 2016; (b) to forecast expected values of aforesaid Islamic 
banking variables, by employing ARIMA models, from January 2017 to December 2019; 
and (c) to investigate the resilience of the Islamic banking industry under different liquidity 
withdrawal scenarios. 

It is pertinent to note that ARIMA modelling approach combines two different processes, 
Autoregressive (AR) and Moving Average (MA), into a single equation. An AR process ex- 
presses a dependent variable Yt as a function of the dependent variable’s past values, while, 
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MA process expresses the dependent variable Yt as a function of the error term’s past values 
(Studentmund, 2016). 
An ARIMA (p, d, g) model can be expressed as follows: 


Y, = Bo + O1%-1 + O2%-2 +++ + OpNi-p + € + PrEr-1 + G2E1-2°++ + bgEr-q 


Where “p” refers to the number of AR terms, “d” denotes the number of times the time 
series is sigerenced before becoming stationary, and “gq” refers to the number of MA terms. 

According to the Box-Jenkin methodology, an ARIMA process consists of following four 
steps: (1) finding out the appropriate values of p, d, and q; (ii) estimating coefficients of AR 
and MA terms of the chosen model; (iii) checking whether the chosen model fits the data; 
and (iv) producing future forecasts on the basis of ARIMA models (Gujarati, 2004). 


ECONOMETRIC ANALYSIS OF THE LIQUIDITY MANAGEMENT ASPECTS 
OF IBF INDUSTRY 


Theoretical Construction of Models 

Model of Banks’ Behaviour in the Competitive Banking Sector 

According to Ismal (2010), a bank’s profit can be calculated by subtracting its asset side’s 
revenues from its liability side’s expenditures. That is: 


mn=r,L+ry—rpD-—-C(D,L) (1) 


In the equation (1), “z’”’ denotes bank’s profit, “rL” denotes interest on loans, “L” denotes 
total outstanding loans, “7” is denoting the money market rate, “rD” denotes interest on 
deposits, “D” is denoting total deposits and “C”’ is the total costs in managing both deposits 
and loans. However, “M” is the bank’s net money market position, which can be formulated 


as: 


M=(1-a)D-L (2) 


“a 29 


In the equation (2), 
reserve requirement. 
When the two equations (1) and (2) are combined together, we get: 


is denoting the reserve coefficient set by a central bank as a cash 


nm(D,L)=(r, -—r)L + [rd - a) -—rp|D -C(D,L) (3) 
A bank’s maximum profit is the first order condition of equation (3). That is: 


On 

OL 
Aforementioned eae (3) and (4) are ee that any competitive banking in- 
stitution will adjust the volume of loans and deposits in such a way that (rL — r) and 
[r(1 — a) — rD] equals its marginal costs of banking operations. Consequently, an increase 
in “rp” will entail a decrease in the bank’s demand for deposits. Similarly, an increase in 


=(7,—-1r)- “(D. L) =0 and = = [rd -a@)-rp]- <(D, L) =0 (4) 
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“r,”’ will entail an increase in the bank’s supply of loans. In case, in an economy, there are 
N number of banks (n = J, 2, ..., N). Each bank is characterised by a loan supply function 
Ln(r,,rp,r) and deposit demand function Dn(rz,rp,r). Let S(rp) is the savings function 
of households and /(r,) is the investment demand by firms. Let banking deposits and T-Bills 
are assumed to be the perfect substitutes for households. The competitive equilibrium of the 
banking sector will be characterised by the following three equations (Freixas & Rochet, 
1999): 


N n 
I(r) = Paes L"(r_.rp,r) (loans market) (5) 


N 
Srp) = B+ ae a D"(rL.rp.1) (savings market) (6) 


N N 
> Eis tp.F) = (alpha) ))_ D"(ri,7p,r) (Interbank market) (7) 


Where ‘B’ is the net supply of T-bills. Equation (7) expresses the fact that aggregate 
position of all banks on the interbank market is zero (M = 0). 

Let us assume that marginal costs of intermediation are constant (C, = yz and Cp = yp) 
then equations (5) and (6) can be replaced by a direct determination of ‘r,’” and ‘rp’, deduced 
from equation (4), such as: r, =r +yz, andrp =r(1-—a@)—- yp. 

Then the interest rate ‘r’ on the interbank market is determined by equation (7), which 
can also be written as: 


L-(@ +41) _ 


S(r(1 = @)yp) - 


B (8) 


e+) =) Luror) =(-@) Der.) (9) 


Above mentioned equations fulfil the expected utility of investors. Equation (8), specifi- 
cally, identifies that liquidity of a bank on its liability side is affected by reserve coefficient 
(qw) and a change in the level of B through open market operations on the equilibrium level 
of interest rates ‘r,’ and ‘rp’. On the other hand, investment demands from firms are in- 
fluenced by money market rate and the costs of intermediation. Equation (9) identifies that 
liquidity of a bank on its asset side is influenced by a set of interest rates (rz, 7p and, r), cost 
of intermediation, total deposits and the reserve coefficient. 


Model of the Liquidity Reserves in Banks 

Liquidity Reserves (R) of a bank include cash and its balances with the central bank as a 
CRR. A bank, therefore, advances credit to its customers equal to the difference between its 
deposits (D) and liquidity reserves (R), that is, DR. Referring to Ismal (2010), let the net 
amount of liquidity withdrawn by a bank is represented by the random variable ‘xX’. When 
the realization ‘x’ of ‘x’ is greater than R, then the bank will experience shortage of liquidity 
and it has to pay a penalty r,(x — R) in proportion of the shortages. 
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In case, banking deposits are assumed as costless and risk free, then the expected profit 
of a bank can be calculated as under: 


mR) = rz,(D — R)+rR-—r,E[max(0, x — R)] (10) 


Let expected costs of liquidity shortages are a convex function and the random variable 
‘x’ has a continuous density function f(x); then these liquidity shortages are differentiable. 
According to Freixas and Rochet (1999), C(R) denotes the expected cost of liquidity re- 
serves, such that: 


C(R) = a) (x — R)f(x)d(x) (11) 

C’(R) = rp | eerie, = -r, Proba [x > R] = 0 (12) 
R 

C”(R) = rp f(R) = 0 (13) 


(Note: C(R) and C(R) are the 1 and 2" difference of C(R)). 
According to Ismal (2010), a bank achieves the maximum profit, when: 


n’(R) = —(r1-—1r)+1rp Proba [X > R] = 0 (14) 


Similarly, the optimum level of liquidity reserves (R*) can be found as follows: 


peas (15) 

lp 
Equation (15) implies that the optimum level of liquidity reserves is that for which the 
marginal opportunity cost of holding reserves becomes equal to the expected cost of liquid- 


ity shortages (Freixas & Rochet, 1999). 


Econometric Analysis 

Liability, Asset, and Liquidity Reserves Models have been estimated econometrically on the 
basis of theoretical construction of these models by employing ARDL modelling approach. 
The aforementioned ARDL models will: (a) explain liquidity related behaviours of Islamic 
banking depositors as well as Islamic bankers; (b) identify those factors that can help in 
determining the asset-liability balance; and (c) identify factors that can help in determining 
the optimal liquidity reserves. 


Definition of Variables and Model Specification 

Liability Model: This ARDL regression model represents the liquidity behaviour of Islamic 
banking depositors. Total Islamic Deposits (TID) is the dependent variable of the model, 
which comprises demand, saving and term deposits of Islamic banking depositors. While, 
independent variables of the model, which are selected on the basis of equation (8), are 
given as under: 
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i. Total Return on Deposits (ROD) for Islamic banking depositors. It helps Islamic bank- 
ing depositors in deciding either to add liquidity or withdraw it form Islamic banks. 

ii. Total Return on Financing (ROF) earned by Islamic banks. It stands as an indicator of 
Islamic banks’ performance and optimal portfolio management. 

iii. Total Cost of Banking Operations (CBO). It stands as a performance indicator of 
robust Islamic bank. 

iv. KIBOR used as a benchmark rate. If KIBOR is higher than ROD paid by Islamic 
banks, return-seeking Islamic banking depositors might transfer their funds from Islamic to 
conventional banks that could cause displaced commercial risk. 

v. The lag of TID. This lag variable is representing the self-assessment of Islamic banking 
depositors regarding their prior investment decisions, and it also influences their liquidity 
behaviour as well. 











TABLE 1 
Descriptive statistics of the liability, asset & liquidity reserves models’ variables 
Variables Mean Median St. Dev Min Max 
TID* 613,290 492,000 456,548 66,011 1,573,000 


Total Islamic Financing (TIF)* 266,861 195,250 195,255 52,621 821,000 
Total Liquidity Reserves (TLR)* 53,100 39,473 35,678 12,165 133,460 


ROD* 17,426 11,220 13,605 2,375 52,392 
ROF* 33,203 22,405 26,108 4,560 99,766 
CBO* 13,244 9,201 10,461 1,805 43,760 
KIBOR** 10.53 10.18 Dene 6.04 15.52 








*Million rupees; ** per cent 


Asset Model 

This ARDL regression model represents the Islamic banks’ liquidity behaviour while deal- 
ing with entrepreneurs on their asset side. TIF is the dependent variable of the model. 
While, independent variables of the model, which are selected on the basis of equation (9), 
are given as under: 

i. ROD received by Islamic banking depositors. This variable is considered by Islamic 
banks while making their portfolio financing decisions in order to retain the existing Islamic 
banking depositors as well as attracting new ones. 

ii. ROF earned by Islamic banks. It stands as an indicator of Islamic banks’ performance 
and a measure of the robustness of its portfolio financing decisions. 

ili. CBO. This variable determines the profitability of Islamic banks. 

iv. KIBOR used as a money market rate. 

v. The lag of TIF. This lag variable is representing the assessment of Islamic banks’ prior 
portfolio financing decisions, and it also influences their liquidity behaviour as well. 


Liquidity Reserves Model 
This ARDL regression model represents the optimal level of liquidity reserves of an Islamic 
bank and those variables that affect the optimum liquidity reserves. TLR is the dependent 
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variable of the model, while, independent variables of the model, which are selected on the 
basis of equation (15), are given as under: 

i. TIF, which influences the optimal level of liquidity reserves of an Islamic bank. 

iit. ROF earned by Islamic banks. This variable is an indicator of Islamic banks robust 
portfolio financing and how attractive and competitive will be the profit sharing to Islamic 
banking depositors. This variable further explains liquidity behaviours of Islamic banks and 
their depositors. 

iii. KIBOR used as a money market rate. This variable influences investment decisions 
of Islamic banking depositors. 

iv. The lag of TLR. This lag variable influences Islamic banks’ reserves. It contains the 
effect of penalty rate (r,) in case of liquidity shortage below the required level of reserves 
and the opportunity cost of holding cash reserves. 


Construction of the Liability, Asset and Liquidity Reserves Models 

Stationarity Test: Since ARDL regression models are estimated using OLS estimation tech- 
niques, therefore, unit root tests are first of all applied to check the stationarity of each 
dependent and independent time series variable. For the purpose, Augmented Dickey and 
Fuller (1979) test is applied. It is interesting to note that Augmented Dickey and Fuller 
(1979) test investigates the unit root, and thus non-stationarity in time series variable. The 
basic idea of the unit root test is explained as follows with the help of a first-order autore- 
gressive model AR (1): 


Y= ut py-1t+& (16) 


Where ‘yw’ and “o? are parameters and ‘e,’ is the white noise. Y, is considered as sta- 
tionary series, if (-—1 < p < 1); similarly, If o = 1, Y; is a non-stationary series. The null 
hypothesis of the Augmented Dickey and Fuller (1979) test is taken as (Ho: p = 1), which 
is tested against the alternate hypothesis (H1: p < 1). 

Augmented Dickey and Fuller (1979) test is carried out by estimating an equation with 
Y,_; subtracted from both sides of the equation (19): 


Y-Yp= et poyi-NYite 


AY, = ut+(p-W)Y-1+ & 


AY, = e+ yY%-1 + & (17) 


Where, AY, is representing changes in the each dependent and independent variable of the 
Asset, Liability, and Liquidity Reserves models; and y = pl. Whereas, corresponding null 
and alternative hypotheses are, Ho: y = 0 and H1: y < 0. 





The parameter ‘p’ is the first-order serial correlation coefficient. 
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The test statistic of the Augmented Dickey and Fuller (1979) test is the f-statistic. The 
null hypothesis of a unit root can be rejected against the alternative, if the t-statistic is less 
than the MacKinnon (1991) critical values for rejection of the aforesaid null hypothesis; 
the time series under examination will then be stationary. Alternatively, the time series will 
be non-stationary, if the f-statistic is greater than the MacKinnon critical values. Akaike 
Information Criterion (AIC) is used for selecting the optimal lag length. 














TABLE 2 
ADF test statistics of liability, asset & liquidity reserves models’ variables 
N=42 Level 1‘ Difference 2" Difference 
t-Stat. Prob.* Test Critical f-Stat. Prob.* Test Critical f-Stat. Prob.* Test Critical 
Values** Values** Values** 
TID 2.322939 1.0000 -4.273277 -2.18822 0.4792 -4.284580 -4.84745 0.0027 -4.296729 
-3.557759 -3.562882 -3.568379 
-3.212361 -3.215267 -3,.218382 
TIF 0.429480 0.9985 -4.262735 0.736817 0.9995 -4.262735 -8.51669 0.0000 -4.262735 
-3.552973 -3.552973 -3,552973 
-3.209642 -3.209642 -3.209642 
TLR -3.56692 0.0454 -4.198503 -7.414310 0.0000 -4.211868 -5.31753 0.0007 -4.252879 
-3.523623 -3.529758 -3,548490 
-3.192902 -3.196411 -3.207094 
ROD -2.20056 0.4753 -4.226815 -3.692845 0.0354 -4.226815 -43.4265 0.0000 -4.226815 
-3.536601 -3.536601 -3.536601 
-3.200320 -3.200320 -3.200320 
ROF -2.08868 0.5349 -4.226815 -3.323780 0.0782 -4.226815 -47.6156 0.0000 -4.226815 
-3.536601 -3.536601 -3.536601 
-3.200320 -3.200320 -3.200320 
CBO -0.17284 0.9914 -4.219126 -3.105372 0.1200 -4.226815 -53.1206 0.0000 -4.226815 
-3.533083 -3.536601 -3.536601 
-3.198312 -3.200320 -3.200320 
KIBOR_ -1.88652 0.6433 -4.198503 -5.389357 0.0004 -4.205004 -5.503944 0.0003 -4.219126 
-3.523623 -3.526609 -3.533083 
-3.192902 -3.194611 -3.198312 








* MacKinnon (1991) one-sided p-values 
**Test Critical Values are basically MacKinnon (1991) critical values for rejection of null hypothesis, which are calculated at 1% level, 
5% level, and 10% level respectively. 


Augmented Dickey and Fuller (1979) test (with an intercept and trend model) statistics re- 
veal dependent variables of Asset, Liability, and Liquidity Reserves models: (1) TID is non- 
stationary at level and 1‘ difference but stationary at 2" difference; (2) TIF is non-stationary 
at level and 1* difference, but stationary at 24 difference; and (3) TLR is non-stationary at 
level (1% statistical significance) but stationary at 1‘ difference. Moreover, independent 
variables of Asset, Liability, and Liquidity Reserves models: (1) ROD is non-stationary 
at level and 1* difference (1% statistical significance) but stationary at 2st difference; (2) 
ROF is non-stationary at level and 1“ difference (1% & 5% statistical significances) but are 
stationary at 2™ difference; (3) CBO is non-stationary at level and 1‘ difference but are 
stationary at 2" difference; and (4) KIBOR is non-stationary at level but stationary at 1‘ 
difference. 
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Since an ARDL model may include I (1) and I (0) variables (but not I (2) variables), 
therefore, the current study integrates all variables in order 1 (1 difference) to find the ro- 
bust Asset, Liability and Liquidity Reserves models. 


Estimation of ARDL Models Using Bounds Testing Approach to Cointegration: This the- 
sis estimate ARDL models by using bounds testing approach to cointegration developed by 
Pesaran, Shin, and Smith (2001). ARDL bounds testing approach involves the following 
two steps (Pesaran & Pesaran, 2009). Firstly, a conditional Error Correction Model (ECM) 
is developed to examine the cointegration among time series variables of interest through an 
F-test, as follows: 


24 
AY, =C+mYitmYit >) YAZi + wAX, +e; (18) 


Where Z,(Y,, X;) Presence of the cointegration between the two time series variables X; and 
Y, can be checked by calculating the F-statistic for the null hypothesis (Ho : 2) = 7, = 0). 
If F-statistic is greater than the upper bound I (1), then null hypothesis HO can be rejected, 
and the alternative hypothesis of long run relationship in levels or cointegration between 
variables is accepted. If, however, F-statistic is less than the lower bound I (0), then null 
hypothesis HO cannot be rejected. In case, the value of F-statistic falls between the lower 
and upper bounds, then results are inconclusive. 

If the two time series variables X; and Y, are cointegrated, then ARDL model can be 
determined with optimal lags selected on the basis of Akaike Information Criterion (AIC) 
in the second step, as follows: 


DL, po)¥Y = Co + B(L, q)X + br (19) 


Where 
O(L, po) = 1— OL + OL? —--» —+@,9L” 


And 
BL, q1) = Bo+ Bil + Pol? +--+ + Ba L® 
Here, L is the lag operator, such that: L'X = X;, —i. 
In the long run, we have ¥, — ¥_-, = +--+ = Y_-po and X; — X;-) = +++ = X;~g1. Thus, we 
can derive the long run cointegration equation as follows: 
Y=anu+tBPxt+y, (20) 
Where 


a = Co/®(1, po); BBC, go); and v, = ,/®C1, po) 


Next, following Banerjee, Dolado, Galbraith, and Hendry (1993) a dynamic ECM for the 
short term relationship can be obtained by linearly transforming the ARDL model. This 
ECM model explains short run dynamics with the long run equilibrium without losing the 
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long run information. Thus, equation (19) can be rearranged as follows to get the short term 
dynamics of the ARDL: 


po-1 gi-1 
AY, =Aao+ >), BpANit )), ByiAXii+ YECTi1 + oh (21) 


Here, the Error Correction Term (ECT;) denotes the OLS residuals series from the long 
run cointegrating regression, which is given as follows: 


ECT, = ¥,- ao - B°X; (22) 


In the equation (21), y the coefficient of ECT;_;, denotes the speed of adjustment, which 
suggests how quickly variables return to equilibrium. Importantly, “y’ should be statistically 
significant and negative (Pak & Ku, 2017). 














TABLE 3 
ADF test statistics of liability, asset & liquidity reserves models’ variables 
ARDL F-Bounds Test t-Bounds Test 
Models 
F-statistic Sig. LowerI(O) UpperI(1) ¢-Statistic Sig. Lower I(0) Upper I(1) 
Liability 10.97449 10% 3.03 4.06 -4.85291 10% -3.13 -4.04 
Model 
5% 3.47 4.57 5% -3.41 -4.36 
2.5% 3.89 5.07 2.5% -3.65 -4.62 
1% 44 372 1% -3.96 -4.96 
Asset 35.80920 10% 3.03 4.06 -2.53123 10% _ -3.13 -4.04 
Model 
5% 3.47 4.57 5% -3.41 -4.36 
2.5% 3.89 5.07 2.5% -3.65 -4.62 
1% 44 S372 1% -3.96 -4.96 
Liquidity 16.51229 10% 3.47 4.45 -8.07883 10%  -3.13 -3.84 
Reserves 
Model 
5% 4.01 5.07 5% -3.41 -4.16 
2.5% 4.52 5.62 2.5% -3.65 -4.42 
1% 5.17 6.36 1% -3.96 -4.73 








Source: I (0) and I (1) asymptotic critical values are obtained from Pesaran et al. (2001) 
Case 5: Unrestricted Constant and Unrestricted Trend. 


Table 3 reveals: (a) In case of liability model, F-bounds test statistic is greater than 
the 1% critical value for the upper bound. Therefore, the hypothesis of “No Long-Run 
Relationship” is strongly rejected. This means that the dependent variable and the four 
regressors are cointegrated. Furthermore, the absolute value of t-bounds test statistic is 
greater than the 2.5% critical value for the upper bound. Therefore, the null hypothesis of 
the t-bounds test is rejected. Hence, the cointegration is not “Nonsensical”. This means 
that the cointegrating relationship is either “usual” or it is “degenerate cointegration”. (b) In 
case of asset model, F-bounds test statistic is greater than the 1% critical value for the upper 
bound. Therefore, the hypothesis of “No Long-Run Relationship” is strongly rejected. This 
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means that the dependent variable and the four regressors are cointegrated. Furthermore, 
the absolute value of t-bounds test statistic is less than the 10% critical value for the lower 
bound. Therefore, the null hypothesis of the t-bounds test cannot be rejected. This means 
that the cointegrating relationship is “Nonsensical”. (c) In case of liquidity reserves model, 
F-bounds test statistic is greater than the 1% critical value for the upper bound. Therefore, 
the hypothesis of “No Long-Run Relationship” is strongly rejected. This shows that the 
dependent variable and the three regressors are cointegrated. Furthermore, the absolute 
value of ¢-bounds test statistic is greater than the 1% critical value for the upper bound. 
Therefore, the null hypothesis of the t-bounds test is rejected. Hence, the cointegration is 
not “Nonsensical”. This means that the cointegrating relationship is either “usual” or it is 
“degenerate cointegration”. 

In conclusion, cointegration between the dependent variables and regressors in the lia- 
bility, asset, and liquidity reserves models have been identified. Therefore, the underlying 
ARDL models can now be established to investigate the long run slope-estimated coeffi- 
cients and the short-run dynamic-estimated coefficients (Pak & Ku, 2017). 











TABLE 4 
ADF test statistics of liability, asset & liquidity reserves models’ variables 
ARDL Models Variable Coefficient Std. Error t-Statistic p-value 
Liability Model D(ROD) -1.235118 8.739515 -0.141326 0.8911 
D(ROF) 0.045116 4.935543 0.009141 0.9929 
D(CBO) 7.136819 9.343495 0.763828 0.4669 
D(KIBOR) 8949.219 7862.807 1.138171 0.2880 
Asset Model D(ROD) -10.37344 13.05325 -0.794702 0.4343 
D(ROF) 1.758696 8.589283 0.204755 0.8394 
D(CBO) 24.62332 11.95876 2.059020 0.0501 
D(KIBOR) -10961.52 7611.455 -1.440134 0.1622 
Liquidity Reserves D(TIF) -0.036908 0.037291 -0.989730 0.3380 
Model 
D(ROF) -1.600123 0.291504 -5.489192 0.0001 
D(KIBOR) 3831.974 1270.014 3.017269 0.0087 








Source: I (0) and I (1) asymptotic critical values are obtained from Pesaran et al. (2001) for 
Case 5: Unrestricted Constant and Unrestricted Trend. 
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Short run error correction estimates of ARDL (3, 4, 5, 5, 5) model 


2018 








Dependent Variable: D(TID) 
(Case 5: Unrestricted Constant and Unrestricted Trend) 








Independent Variable Coefficient f-statistic p-value 
Constant -4975.916 -1.255131 0.2449 
Trend 2438.827 7.661414 0.0001 
D(TID(-1)) 0.657200 4.329086 0.0025 
D(TID(-2)) 0.969765 5.061428 0.0010 
D(ROD) 1.996523 0.993472 0.3496 
D(ROD(-1)) -8.638917 -3.233833 0.0120 
D(ROD(-2)) -12.78015 -3.249548 0.0117 
D(ROD(-3)) -7.475340  -2.106261 0.0683 
D(ROF) -1.120193  -0.985833 0.3531 
D(ROF(-1)) 5.018159 2.995356 0.0172 
D(ROF(-2)) 5.950740 2.340946 0.0473 
D(ROF(-3)) 2.125687 0.960049 0.3651 
D(ROF(-4)) -1.970427 -4.190991 0.0030 
D(CBO) 0.080961 0.036407 0.9718 
D(CBO(-1)) -7.137212  -2.237163 0.0557 
D(CBO(-2)) -1.526826 -0.448123 0.6660 
D(CBO(-3)) 3.598946 1.268418 0.2403 
D(CBO(-4)) 6.130394 3.518699 0.0079 
D(KIBOR) -85.32015  -0.040346 0.9688 
D(KIBOR(-1)) -12975.22  -5.115744 0.0009 
D(KIBOR(-2)) -11252.34 -4.185489 0.0031 
D(KIBOR(-3)) -13539.56 -4.226890 0.0029 
D(KIBOR(-4)) -6900.661 -2.925284 0.0191 
ECTt-1 -1.506994 -9.072412 0.0000 
Diagnostic Analysis Value p-value 

R-squared 0.995245 

Residual Sum of Square 1.00E+09 

AIC DA Wed bee!) 

F-Statistics 109.2109 0.0000 
Jarque-Bera Test 14.10160 0.0009 
Breusch-Godfrey Serial Correlation LM Test 0.067103 0.9358 
Heteroscedasticity Test: Breusch-Pagan-Godfrey 0.236394 0.9978 

Ramsey RESET Test 3.260103 =0.1140 
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TABLE 6 


Short run error correction estimates of ARDL (1, 2, 2, 2, 1) model 








Dependent Variable: D(TIF) 
(Case 5: Unrestricted Constant and Unrestricted Trend) 








Independent Variable Coefficient f-statistic p-value 
Constant -1032.693 -0.281749 0.7805 
Trend 188.6381 1.243207 0.2253 
D(ROD) 0.092500 0.058182 0.9541 
D(ROD(-1)) 8.095180 5.004351 0.0000 
D(ROF) -0.732756 -0.792951 0.4353 
D(ROF(-1)) -3.918659 -4.290403 0.0002 
D(CBO) 4.995538 6.790088 0.0000 
D(CBO(-1)) -3.704249 -5.264724 0.0000 
D(KIBOR) -2559.633 -1.789674 0.0856 
ECTt-1 -0.457431 -14.41157 0.0000 
Diagnostic Analysis Value p-value 

R-squared 0.953782 

Residual Sum of Square 3.05E+09 

AIC 21.52415 

F-Statistics 66.49630 0.0000 
Jarque-Bera Test 1.643627 0.4396 
Breusch-Godfrey Serial 0.695579 0.5090 

Correlation LM Test 

Heteroscedasticity Test: 1.917215 0.0787 
Breusch-Pagan-Godfrey 

Ramsey RESET Test 5.809881 0.0240 
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TABLE 7 


Short run error correction estimates of ARDL (4, 4, 3, 5) model 


2018 








Dependent Variable: D(TLR) 


(Case 5: Unrestricted Constant and Unrestricted Trend) 








Independent Variable Coefficient f-statistic p-value 
Constant -12559.48 = -2.958447 0.0098 
Trend 1906.409 7.148486 0.0000 
D(TLR(-1)) 3.175308 6.443039 0.0000 
D(TLR (-2)) 1.985704 5.709212 0.0000 
D(TLR (-3)) 0.966628 5.392487 0.0001 
D(TIF) -0.339112 -2.524189 0.0234 
D(TIF(-1)) -0.021172 -0.127632 0.9001 
D(TIF(-2)) 0.635507 3.877317 0.0015 
D(TIF(-3)) 0.396665 2.186007 0.0451 
D(ROF) -1.789084 -7.914087 0.0000 
D(ROF(-1)) 4.812320 9.019295 0.0000 
D(ROF(-2)) 2.500677 8.646634 0.0000 
D(KIBOR) 5409.922 3.206298 0.0059 
D(KIBOR(-1)) -10802.37  -5.293423 0.0001 
D(KIBOR(-2)) -4956.414 -2.784185 0.0139 
D(KIBOR(-3)) -9229.545 -4.683715 0.0003 
D(KIBOR(-4)) -6212.491 -3.543796 0.0029 
ECT,- -5.337371 -8.902752 0.0000 
Diagnostic Analysis Value p-value 

R-squared 0.937858 

Residual Sum of Square 1.64E+09 

AIC 21.47066 

F-Statistics 15.98005 0.0000 
Jarque-Bera Test 9.961281 0.0069 
Breusch-Godfrey Serial Correlation LM Test 0.291616 0.7518 
Heteroscedasticity Test: Breusch-Pagan-Godfrey 0.819465 0.6668 

Ramsey RESET Test 1.433827 0.2510 








Diagnostic Analysis: In order to assess the validity of the assumptions that underlie the li- 
ability, asset, and liquidity reserves models, different diagnostic tests have been conducted, 
whose results are shown in Tables 5, 6, and 7 respectively. Breusch-Godfrey serial corre- 
lation LM test statistic confirms that residuals in the liability, asset and liquidity reserves 
models are serially uncorrelated. Moreover, Breusch-Pagan-Godfrey heteroscedasticity test 
statistic reveals that residuals in the aforementioned ARDL models are homoscedastic. In- 
terestingly, functional form of the aforesaid three ARDL models is checked by applying 
Ramsey RESET test. Results of the said test confirm that there is no misspecification with 
respect to the functional form of the three ARDL models. In order to examine the stability 
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of the estimated parameters of the said ARDL models, CUSUM tests are applied. Figures 1, 
2, and 3 respectively show that all the three estimated ARDL models are valid with robust 
recursive residuals. 
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FIGURE 1. CUSUM of recursive residuals of liability model 
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FIGURE 2. CUSUM of recursive residuals of asset model 


72 Rafiullah, S., Khan, A. & Mumtaz, F: R. - Liquidity management by Islamic .... 2018 


12 


2013 2014 2015 2016 


——  CUSUM ---— 5% Significance 


FIGURE 3. CUSUM of recursive residuals of liquidity reserves model 


Interpretations of the Models 

Liability Model 

Islamic banks’ liability side has a specific liquidity behaviour that represents Islamic bank- 
ing depositors’ liquidity behaviour. Firstly, the most important factor that determines the 
liquidity behaviour of Islamic banking depositors is the amount of returns paid to them in 
a short-term period. It is interesting to note that higher the Islamic banks’ profit (return on 
financing), better the Islamic banking depositors’ perception of their banks that will ulti- 
mately lead them to deposit more funds. 

Secondly, after returns on deposits, the amount and level of benefits received by Islamic 
banking depositors by depositing money stimulate their investment decision. Besides mon- 
etary benefits (returns, bonuses, etc.), other non-monetary benefits, for instance, standard 
banking facilities, use of innovative technologies, branch network, etc., improve Islamic 
banking depositors satisfaction and loyalty. Interestingly, more satisfied and loyal Islamic 
banking depositors will add more funds in their banking deposits. 

Thirdly, KIBOR-a benchmark rate-is an important factor explaining Islamic banking de- 
positors’ behaviour. Therefore, for being competitive so as to convince depositors to add 
more funds, returns on deposits must be competitive and matching the expected return of 
Islamic banking depositors (KIBOR). 

Fourthly, costs of banking operations also predict Islamic banking depositors’ behaviour. 
Higher the costs of banking operations, lower will be the cost efficiency and lesser will be 
the amount available for distributing among Islamic banking depositors. This will provoke 
Islamic banking depositors not to add funds into their existing deposits, may be causing 
withdrawal of funds. 


Asset Model 

The asset side of Islamic banks also has a specific liquidity behaviour explaining the liq- 
uidity behaviour of Islamic banks. Firstly, one of the Islamic banks’ purposes in managing 
liquidity is optimizing returns on financing. If such returns tend to be more prospective, 
Islamic banks will advance more financing. 
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Secondly, returns on deposits paid to Islamic banking depositors is the factor affecting 
Islamic financing. This factor is important because it is related to the Islamic banks’ com- 
petitiveness, attractiveness and their ability to match depositors’ expectations. In this way, 
Islamic banks will attract more deposits and, therefore, more funds will be available to fi- 
nance the projects. 

The third variable that affects Islamic financing is the KIBOR, as a benchmark rate. If 
such rate is high, Islamic banks tend to advance more financing. Whereas, Islamic banks 
total financing in the previous periods is the fourth variable that influences their liquidity 
behaviour, that also guides their future financing policy. This variable estimates the success 
or failure of Islamic banks’ portfolio financing policy and suggests future policy regarding 
portfolio financing. 

Finally, the fifth variable affecting Islamic banks’ financing is the cost of banking oper- 
ations, which is associated with the Islamic banks’ efforts in evaluating the performance of 
their business partners and in finding the prospective projects to be financed. Sustainable 
and higher profits can be made by: (a) selecting financially stable entrepreneurs as business 
partners; (b) analysing critically their business proposals, and (c) building amicable working 
relations with business partners. 

When Islamic banks expand financing, associated cost might also increase accordingly. 
However, if returns on financing are as expected, Islamic banks may achieve more returns by 
expanding financing and increasing associated costs. Indeed, these additional costs of bank- 
ing operations will benefit Islamic banks in competing conventional banks in a dual banking 
system, attracting more deposits, and fulfilling expectations of existing Islamic banking de- 
positors. 


Liquidity Reserves Model 
This ARDL regression model discovers the factors, as indicated here, that determine the 
optimal level of liquidity reserves. Firstly, the current liquidity reserves position depends on 
returns on financing. If returns on financing fall then the liquidity reserves model suggest 
that Islamic banks will need more liquidity since the returns on deposits paid to Islamic 
banking depositors might as well fall and, therefore, return-seeking depositors might with- 
draw their deposits. Secondly, liquidity reserves depend on KIBOR. If the return paid to 
Islamic banking depositors fall short of KIBOR, then Islamic banks will need more liquid 
funds in anticipation of displaced commercial risk. 

The third factor that affects the liquidity reserves position is the Islamic banks’ financing. 
If it goes up, Islamic banks need to raise the liquidity reserves level. It is interesting to note 
that expansion in financing is basically possible when Islamic banking depositors add more 
funds. These high deposits require Islamic banks to hold higher liquidity reserves (Cash 
Reserve Requirement of the central bank). Finally, the previous liquidity reserves’ position 
explains the pattern of liquidity demands by Islamic banking depositors and entrepreneurs. 


Resilience of the Industry against Liquidity Pressures 
The current study tries to investigate the resilience of the Islamic banking industry by 
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analysing: (a) the ability of liquidity suppliers to meet liquidity demands; and (b) the level 
of Islamic banking depositors’ liquidity withdrawals, which could possibly make an Islamic 
bank fail to meet its financial obligations. For the purpose, ARIMA modelling approach is 
employed: (a) to estimate models for variables related to Islamic bank’s liquidity manage- 
ment by employing historical data from July 2006 to December 2016; (b) to forecast ex- 
pected values of aforesaid Islamic banking variables, by employing ARIMA models, from 
January 2017 to December 2019; and (c) to investigate the Islamic banking industry’s re- 
silience under different liquidity withdrawal scenarios. 


ARIMA Process 

ARIMA modelling process involves eight variables, which are sorted into following two 
groups: (1) Liquidity-demanders-Demand Deposits (DD), Saving Deposits (SD), Term De- 
posits (TD); and (2) Liquidity-suppliers TLR that include cash and balances with treasury 
banks, Interbank Placements (IP), Shari‘ah-Compliant Investments (SI), Equity Participa- 
tion (EP), and Funds Borrowed (FB) from IFIs. Variables of liquidity-supplier group will 
be subdivided into: (a) 1*' Tier Liquid Instrument-that serve the liquidity demands from 
Demand Deposits and Saving Deposits; and (b) 2" Tier Liquid Instrument-that serve all 
liquidity demands from Demand, Saving and Term Deposits. The current study considers 
TLR and FB as 1“ tier liquid instruments, while, TLR, FB, IP, SI, and EP collectively are 
considered as 2" tier liquid instruments. 


Identification of Variables: Firstly, descriptive statistics of liquidity-demanders and liquidity- 
suppliers group of variables are given in the table 8. Next, each variable of liquidity- 
demanders (DD, SD, TD) and liquidity-suppliers (TLR, IP, SI, EP, and FB) are tested for 
stationarity. Table 9 provides the Augmented Dickey and Fuller (1979) test statistics of each 
variable of liquidity-demanders and -suppliers. 


TABLE 8 
Descriptive statistics of liquidity-demanders & liquidity-suppliers 








Liquidity-Demanders 








Variables Mean Median St. Dev Min Max 
DD* 157,857 110,356 137,573 13,379 483,434 
SD* 222,512 160,302 183,844 18,684 633,664 
TD* 179,410 175,762 106,277 24,399 362,530 
Liquidity-Suppliers 

Variables Mean Median St. Dev Min Max 
TLR* 58,410 43,420 39,246 13,382 146,806 
Interbank Placements (IP)* 27,558 26,042 10,847 13,113 61,362 
SI* 234,495 228,323 192,252 3,437 661,787 
EP* 52,626 53,029 26,380 13,542 112,223 
FB from IFIs * 189,387 143,356 137,885 36,456 579,128 








* Million rupees 
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TABLE 9 
ADF test statistics of liquidity-demanders and liquiditysSuppliers 


75 








Liquidity-Demanders 











N=42 Level 1° Difference 2™ Difference 
t-Stat. Prob.* Test Critical t-Stat. Prob.* Test Critical t-Stat. Prob.* Test Critical 

Values** Values** Values** 

DD -0.44445 0.9819 -4.226815 -1.84723 0.6611 -4.226815 -8.53478 0.0000 -4.226815 
-3.536601 -3.536601 -3.536601 
-3.200320 -3.200320 -3.200320 

SD -0.02391 0.9946 -4.198503 -7.19718 0.0000 -4.205004 -8.06577 0.0000 -4.219126 
-3.523623 -3.526609 -3.533083 
-3.192902 -3.194611 -3.198312 

TD -3.74343 0.0310 -4.211868 -4.56087 0.0047 -4.252879 -4.31058 0.0097 -4.296729 
-3.529758 -3.548490 -3.568379 
-3.196411 -3.207094 -3.218382 

Liquidity-Suppliers 
N=42 Level 1°‘ Difference 2™ Difference 
t-Stat. Prob.* Test Critical t-Stat. Prob.* Test Critical t-Stat. Prob.* Test Critical 

Values** Values** Values** 

TLR 4.283972 1.0000 -4.198503 -2.83658 0.1934 -4.205004 -8.67464 0.0000 -4.211868 
-3.523623 -3.526609 -3.529758 
-3.192902 -3.194611 -3.196411 

IP -1.71747 0.7253 -4.198503 -5.23074 0.0006 -4.205004 -9.36107 0.0000 -4.211868 
-3.523623 -3.526609 -3.529758 
-3.192902 -3.194611 -3.196411 

SI -2.50215 0.3255 -4.205004 -3.48416 0.0548 -4.205004 -7.36613 0.0000 -4.211868 
-3.526609 -3.526609 -3.529758 
-3.194611 -3.194611 -3.196411 

EP -4.46470 0.0059 -4.252879 -4.93262 0.0014 -4.205004 -7.23344 0.0000 -4.219126 
-3.548490 -3.526609 -3.533083 
-3.207094 -3.194611 -3.198312 

FB -0.94640 0.9404 -4.198503 -4.75274 0.0024 -4.211868 -6.82530 0.0000 -4.219126 
-3.523623 -3.529758 -3.533083 
-3.192902 -3.196411 -3.198312 








* MacKinnon (1991) one-sided p-values. 
**Test Critical Values are basically MacKinnon (1991) critical values for rejection of null hypothesis, which are calculated at 1% level, 
5% level, and 10% level respectively. 


Table 8 reveals that all variables under investigation are non-stationary at level but sta- 
tionary at 1°‘ difference, except DD, TLR and SI that are stationary at 2" difference. 

The next identification process in the ARIMA modelling is to examine the Autoregres- 
sive and Moving Average patterns using a correlogram test for behaviour patterns of Au- 
tocorrelation Function (ACF) and Partial Correlation Function (PACF) patterns. The three 


commonly found patterns in the ARIMA models are as follows: 


a. All ACF values are zero; 
b. ACF values show a cut off pattern between 1‘ and the 2™ or 3" periods of ACF; and 


c. ACF values are showing a dying down pattern. 


In ARIMA modelling, if ACF is dying down and PACE cuts off, then following pure AR 
model should be employed: 
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Li = re) + 6, Z;-1 + 6523-2 + +6,Z;—p + Et (23) 


In case, ACF cuts off and PACF is dying down, then following pure Moving Average 
(MA) model should be employed: 


Z: = e+ & — Diet —1- Qt —2—----— Det -—q (24) 


In case, ACF and PACF both are dying down, then following model, which is a combina- 
tion of AR and MA models, should be employed: 


Zr = U+ Oj Zy-1 + O2Z;-2 + +++ + OyLZy_p + &; — DyEt — 1 — O9eEt —2—---— Det — g (25) 


Table 10 reveals that: (a) for DD, ACF is dying down, while, PACF cuts off; (b) for SD, 
ACF and PACF both are dying down; (c) for TD, ACF cuts off, while, PACF is dying down; 
(d) for TLR, ACF cuts off, while, PACF is dying down; (e) for IP, ACF cuts off, while, PACF 
is dying down; (f) for SI, ACF is dying down, while, PACF cuts off; (g) for EP, ACF and 
PACF both are dying down; and (h) for FB, ACF and PACF both are dying down. 








TABLE 10 
ACF & PACF correlograms 
DD SD TD TLR IP SI EP FB 





Lags ACF PACF ACF PACF ACF PACF ACF PACF ACF PACF ACF PACF ACF PACF ACF PACF 





1 0.91 0.91 0.92 0.92 0.94 0.94 0.90 0.90 0.84 0.84 0.91 0.91 0.90 0.90 0.89 0.89 
2 0.82 -0.05 0.84 -0.01 0.88 -0.07 0.80 -0.03 0.62 -0.26 0.81 -0.10 0.78 -0.11 0.76 -0.20 
3 0.74 -0.02 0.77 -0.02 0.81 -0.06 0.70 -0.03 0.37 -0.22 0.71 -0.01 0.66 -0.10 0.64 0.07 
4 0.66 -0.05 0.68 -0.10 0.74 -0.12 0.61 -0.04 0.22 0.20 0.65 0.12 0.53 -0.08 0.56 0.03 
5 0.58 -0.03 0.61 -0.01 0.66 -0.02 0.53 -0.01 0.18 0.22 0.59 -0.03 0.42 -0.02 0.48 -0.05 
6 0.50 -0.04 0.53 -0.03 0.58 -0.10 045 -0.02 0.18 -0.11 0.54 -0.02 0.32 -0.00 0.42 0.06 
7 0.43 -0.01 0.46 0.00 O51 0.01 0.38 -0.01 0.20 -0.01 0.48 0.00 0.24 -0.01 0.37 -0.00 
8 0.36 -0.04 0.40 -0.06 044 0.05 0.32 -0.02 0.21 0.13 043 -0.02 0.18 0.00 0.31 -0.05 
9 0.30 -0.03 0.34 0.01 0.38 -0.01 0.26 -0.02 0.18 -0.11 0.38 -0.01 0.16 0.17 0.28 0.10 
10 0.23 -0.04 0.28 -0.02 0.33 -0.02 0.21 -0.02 O11 -0.11 034 -0.01 0.14 -0.03 0.25 -0.03 
11 0.18 0.00 0.22 -0.09 0.27 -0.12 0.16 -0.01 0.07 0.16 0.29 -0.03 0.13 -0.04 0.22 -0.02 
12 O12 -0.04 0.16 -0.03 0.20 -0.13 0.12 0.01 0.02 -0.04 0.24 -0.10 0.12 -0.00 0.18 -0.06 
13 0.07 -0.02 0.10 -0.04 0.14 -0.01 0.09 -0.02 0.03 -0.03 0.17 -0O.11 O11 O01 0.13 -0.03 
14 0.02 -0.03 0.04 -0.06 0.07 -0.05 0.05 -0.03 0.01 -0.12 0.09 -0.06 0.10 -0.01 0.08 -0.06 


15 -0.02 -0.01 -0.01 -0.03 0.01 -0.01 0.02 -0.03 -0.01 0.06 0.02 -0.10 0.08 -0.03 0.03 -0.04 
* MacKinnon (1991) one-sided p-values. 
**Test Critical Values are basically MacKinnon (1991) critical values for rejection of null hypothesis, which are calculated at 1% level, 
5% level, and 10% level respectively. 








Estimation of Models: All the eight estimated ARIMA models are presented below along 
with the values of coefficients, f-statistics (in brackets), R-squared, AIC and F-statistics with 
the p-value in brackets. 


588.5730  -0.799092 = -0.459417 ~— -0.598620 
[3.193323] [-6.496046] [-2.867596] [-4.927881] 
R-squared 0.778280 AIC 19.3099 1 F-stat 30.71422 [0.0000] 
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ASD, = L + OSD;_1 + Et = DM E;-] = MD €;_2 (27) 
778.1852 0.561544 = -1.999983 0.999983 
[15.12540] [3.589414] [-903.6047] [644.4647] 
R-squared 0.670432 AIC 21.95819  F-stat 17.79988 [0.0000] 
AT D, = fot &; — D1 E;-1 — O2E+-2 — D3E;-3 — D4e;-4 (28) 
8324.090 0.254683 -0.275776 = -0.344774 = -0.634128 
[18.66156] [7.310027] [-7.525896] [-7.482580] [-5.898475] 
R-squared 0.368723 AIC 21.62881 F-stat 4.088637 [0.0049] 
161.0087 -0.421087 
[1.147623] [-2.793731] 
R-squared 0.143750 AIC 17.62521 F-stat 3.105842 [0.0566] 
AIP, = + & — D1 &;-1 — D2E;_-2 — O3E;-3 — Dye;_4 (30) 
653.0106 0.070301 0.215291 -0.722382  -0.563196 
[2.329949] [0.000366] [0.001449] [-0.006175] [-0.002402] 
R-squared 0.349521 AIC 19.64460 F-stat 3.761294 [0.0078] 
ASI, = HO, STy-1 + €; (31) 
16192.19 0.538062 
[2.068020] [4.192150] 
R-squared 0.301028 AIC 23.14805 F-stat 8.182768 [0.0011] 
AEP, = L + O, EP, + O,E P;_> + Et a M1 E;-] a DxE;_2 (32) 
2191.220 1.669098  -0.877750 ~— -1.640258 0.999993 
[2.369016] [13.85902] [-8.314473] [-12.23886] [6.947091] 
R-squared 0.297128 AIC 19.53089 F-stat 2.959141 [0.0249] 
AFB, = w+ 0, FB, + €; — Dy E;-1 — O2E;-2 — D3E;-3 — D4e;-4 (33) 
8955.475 = -0.985477 1.556908 0.367426 = -0.036826 0.193013 
[2.339270] [-11.70936] [0.007043] [0.005455] [-0.001147] [0.000966] 


R-squared 0.341640 AIC 21.88614  F-stat 2.940583 [0.0202] 


Forecasting on the Basis of Estimated ARIMA Models: In order to forecast the future liquid- 
ity demands and its supply and assess the Islamic banking industry’s resilience against an 
unanticipated liquidity demands, the eight estimated ARIMA models are utilized to generate 
estimated values from January 2017 to December 2019. 

The current study anticipates the first scenario of regular liquidity withdrawals by Islamic 
banking depositors. The second scenario is of irregular liquidity withdrawals-liquidity de- 
mands rise than the former scenario. This scenario can arise when Islamic banking deposi- 
tors withdraw their deposits due to unfavourable economic conditions and hold more cash. 
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The last scenario is that of liquidity run, which might possibly be occurring due to any 
financial sector’s crisis or due to certain circumstances when Islamic banking depositors’ 
confidence over their banks is shattered. 


Resilience of the 1°t Tier Liquid Instruments 

In order to investigate the 1‘ tier liquid instruments’ resilience, the aforementioned three 
liquidity withdrawal scenarios are examined. The first scenario is assumed to be of regu- 
lar liquidity withdrawals by Islamic banking depositors. The current study assumes regular 
liquidity withdrawals equal to the 10% of each quarterly balance of demand and saving de- 
posits. 1‘ tier liquid instruments’ resilience in case of regular liquidity withdrawals is drawn 
as a thick line in the Figure 4. 

The second scenario is considered as of irregular liquidity withdrawals, which is assumed 
to be equal to 25% of each quarterly balance of demand and saving deposits. 1*' tier liquid 
instruments’ resilience in case of any irregular liquidity withdrawals is drawn as a dotted 
line in the Figure 4. 














Gg Tier 1 Liquid Instruments 

Regular Liquidity Withdrawals 
—+— Irregular Liquidity Withdrawals 
—— Liquidity Run 
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FIGURE 4. Tier-1 liquid instruments’ resilience 


The last scenario of liquidity run is assumed with 45% of each quarterly balance of de- 
mand and saving deposits. A more severe scenario (more than 45%) of liquidity run is not 
taken into consideration because 45% liquidity withdrawals would be giving strong signals 
to Islamic banks to take necessary actions in order to avoid further worsening of the sce- 
nario. 1*' tier liquid instruments’ resilience in case of a liquidity run is drawn as a starred 
line in the Figure 4. 


Resilience of the 2" Tier Liquid Instruments 
24 tier liquid instruments improve the liquidity supply for tackling the liquidity demands 
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from term deposits in addition to the liquidity demands from demand and saving deposits. 
The first scenario is of regular liquidity withdrawals. The current study assumes that 10% 
term deposits of each quarterly balance are terminated in addition to the regular liquidity 
withdrawals equal to the 10% of each quarterly balance of demand and saving deposits. 2™ 
tier liquid instruments’ resilience in case of regular liquidity withdrawals is drawn as a thick 
line in the Figure 5. 
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FIGURE 5. Tier-2 liquid instruments’ resilience 


The second scenario is of the irregular liquidity withdrawals that are assumed to be the 
termination of 25% term deposits of each quarterly balance in addition to the liquidity with- 
drawals from demand and saving deposits equal to the 25% of each quarterly balance. 2" 
tier liquid instruments’ resilience in case of any irregular liquidity withdrawals is drawn as 
a dotted line in the Figure 5. 

The last scenario of liquidity run is assumed with termination of 45% term deposits of 
each quarterly balance in addition to the liquidity withdrawals equal to the 45% of each 
quarterly balance of demand and saving deposits. 2™ tier liquid instruments’ resilience in 
case of liquidity run is drawn as a starred line in the Figure 5. 


Findings of the Resilience Analysis 
Islamic banking industry’s resilience analysis lads to the the following findings that are im- 
portant with regard to the Islamic banks’ liquidity management: 

a) Islamic banking industry performed historically well in managing the liquidity. 

b) Tier-2 liquid instruments performed well historically in mitigating liquidity run condi- 
tions; moreover, it is forecasted that tier-2 liquid instruments would possibly be performing 
well in mitigating any future liquidity run conditions (see Figure 6). Likewise, tier-1 liq- 
uid instruments performed historically well in mitigating liquidity run conditions, however, 
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tier-1 liquid instruments would apparently be facing a liquidity mismatch in the last quarter 
of 2019 (see Figure 7). 
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FIGURE 6. Performance of tier-2 liquid instruments 


c) Based on the severe liquidity run scenarios, tier-1 liquid instruments have absorbed the 
liquidity pressure of up to 55% of withdrawals till 2013. Afterwards, 1*' tier fails to handle 
liquidity run. However, tier-2 liquid instrument have managed to handle liquidity pressure 


of up to 95% of withdrawals till 2013 (see Figure 8). 
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FIGURE 7. Performance of tier-1 liquid instruments 
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FIGURE 8. Performance of liquid instruments in severe liquidity run scenarios 


d) The percentage assumptions of liquidity run delivers the important message that Is- 
lamic banks’ failure in liquidity management begins with this percentage of liquidity with- 
drawals. Following are some of the steps that can be taken to maintain the sound conditions 
of Islamic banking industry and prevent any liquidity runs: (1) educating the Islamic banking 
depositors and general public about the contract mechanics of Islamic banking transactions; 
(ii) improving the asset-liability management; and (111) optimizing portfolio financing in or- 
der to realize better returns on financing and paying competitive returns to Islamic banking 
depositors. 

e) There is another tier (tier-3) of liquid instruments available to Islamic banks in case of 
any emergency. Tier-3 liquid instruments mainly include liquidity provided by State Bank 
of Pakistan as a Lender of Last Resort. However, using tier-3 liquid instruments may have 
certain negative implications for Islamic banks. 


CONCLUSION 


This empirical study was designed to analyse different liquidity management aspects of 
Islamic banking industry of Pakistan. For this purpose, econometric models identifying 
factors that influence Islamic banks in balancing liquidity on liability as well as asset sides 
and in maintaining the optimum liquidity reserves have been developed. Furthermore, the 
study has also endeavoured to investigate the resilience of the Islamic banking industry 
against anticipated and unanticipated liquidity pressures. 

It is pertinent to note that Islamic banking liability and asset models, with regard to liq- 
uidity management, have identified the significant role of following variables: (a) returns 
on deposits, (b) returns on financing, (c) costs of banking operations, and (d) KIBOR. Is- 
lamic banking liquidity reserves model recommends Islamic banks to consider following 
variables, while developing optimum liquidity reserves: (a) total Islamic financing, (b) re- 
turns on financing, and (c) KIBOR. 

Resilience analysis of Islamic banking industry has found that liquid instruments per- 
formed well historically in mitigating liquidity run conditions. Furthermore, it is forecasted 
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that tier-2 liquid instruments would possibly be performing well in mitigating any future liq- 
uidity run conditions (up to 95% of deposits). However, in case of tier-1 liquid instruments, 
there is a possibility of liquidity mismatches when liquidity withdrawals exceed the limit of 
55% of deposits. 

Findings of the current study suggest that Islamic banking depositors, besides their reli- 
gious motives of supporting Islamic banks, expect from their banks to produce profits and 
pay competitive returns on their deposits. Therefore, Islamic banks need to make prudent 
portfolio financing so as to pay competitive returns to their depositors. 

It is pertinent to note that major banking business is that of maturity transformation of 
short-term deposits into long-term project financing. Inability of Islamic banks to meet 
depositors’ liquidity withdrawal demands exposes them to liquidity risk. Consequently, liq- 
uidity management becomes a critical issue that calls for prompt responsiveness by Islamic 
banks and prudential supervision by banking regulators. In this regard, econometric models 
developed in the current study will help Islamic banks in managing deposits, financing port- 
folios, and liquidity reserves, and enable them for forecasting the future liquidity demands 
and supply. 
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